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Abstract 

The results of four recently introduced beamforming schemes for phased array systems are discussed, each 
of which is capable to provide high sensitivity and accurate polarimetric performance of array-based radio 
telescopes. Ideally, a radio polarimeter should recover the actual polarization state of the celestial source, 
and thus compensate for unwanted polarization degradation effects which are intrinsic to the instrument. In 
this paper, we compare the proposed beamforming schemes through an example of a practical phased array 
system (APERTIF prototype) and demonstrate that the optimal beamformer, the max-SLNR beamformer, 
the eigenvector beamformer, and the bi-scalar beamformer are sensitivity equivalent but lead to different 
polarization state solutions, some of which are sub-optimal. 



1 Introduction 

The radio astronomy community has recently decided upon the science cases that future radio telescopes 
must enable, which directly translate into stringent requirements on the instrument's survey speed. These 
include beam sensitivity, size of the field-of-view (FoV), number of simultaneous beams, and frequency 
bandwidth, as well as on the polarization purity of the instrument. Conventional single-beam antennas are 
not capable to conform to these multi-beaming performance specifications. Research has therefore been 
focused on phased- array antenna concepts, since then the individual antenna array output signals can be 
weighted to realize multiple desired beams. Furthermore, array signal processing techniques can be used 
to exploit the degrees of freedom in the beamformer weights for increased beam sidelobe level control, 
interference rejection, and gain and noise level control of the beams (SNR) [HE]. 

To date, the polarimetric response of the antenna array has not received much of attention and only one 
polarization is usually considered at a time. Further, as the radio polarimeter is essentially summing up vector 
element patterns, the polarization characteristics of the realized beam is weight dependent too. In turn, this 
enables us to create two beams in one direction forming an orthogonal polarimetric beam pair. Moreover, 
signal processing techniques can be developed to collectively optimize for the SNR and orthogonality of a 
polarimetric beam pair. This paper assesses the results for four different beamforming schemes that have 
recently been proposed: the optimal beamformer, the max-SLNR beamformer, the eigenvector beamformer, 
and the bi-scalar beamformer. The beam pair sensitivities are shown for the APERFIF prototype system [2], 
for all beamforming scenarios. The polarimetric quality of each pair of nominally orthogonally polarized 
beams is expressed in terms of a Jones matrix. This formulation is used in a numerical model to study the 
performance of several possible algorithms for polarimetric array beamformer weight design. 



2 Polarimetric Beamforming 



Figure [T] provides a system level view of a polarimetric antenna system. The polarimetric properties of the 
source are defined by the 2x2 source covariance matrix Rg, which gives the covariance of the source signal in 
two polarizations defined in a suitable coordinate system. The transfer function of our instrument, including 
antennas, receiver chains and a beamformer, can be described in the form of a 2 x 2 Jones matrix J, which 
describes how the input voltages in two polarizations are transformed into two polarimetric output signals. If 
we correlate these output signals, we get the beamformer output covariance matrix Rq. This is used as input 
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Figure 1: A system level view of a radio polar imeter. 

for further processing, that ultimately produces a reconstructed source covariance matrix Rg that should 
ideally be proportional to Rg. An ideal system does not require polarimetric correction, since it has J = I, 
where I denotes the identity matrix, i.e., it leaves the covariance matrix of the input signal unchanged and 
does not introduce so-called instrumental polarization. By choosing appropriate beamformer weights, the 
instrumental polarization introduced by the antennas and other analog electronics can be compensated for 
in the beamformer. In this paper, we compare the polarimetric and sensitivity performance of four proposed 
beamformers: 

optimal beamformer In this beamforming scheme, we assume that we know the response of the system to 
two reference signals with mutually perfectly orthogonal polarization. We then maximize the sensitivity 
under the constraint that J = I by minimizing the receiver output noise power [3^. 

max-SLNR beamformer The optimal beamformer imposes a hard constraint on the instrumental transfer 
function. By optimizing the signal-to-leakage-and-noise ratio, one can turn this into a soft constraint 
that does not suppress the leakage further than the noise level in the calibration measurement. This 
ensures that the noise in the calibration measurement is not over-interpreted. 

eigenvector beamformer This more practical beamformer is based on the fact that an unpolarized signal 
can be regarded as a superposition of two uncorrelated orthogonally polarized signals. We can thus 
produce two orthogonal output signals by measuring the signal covariance matrix of the antenna 
outputs, which is a rank-2 matrix, and associating each beamformer output with one of the dominant 
eigenvectors [4]. If the phased array feed consists of two sets of single polarization elements, we can 
make an approximate polarimetric correction by exploiting the intrinsic polarization of each set of 
elements. 

bi-scalar beamformer This beamformer assumes that the array consists of two sets of single polariza- 
tion elements and simply applies max-SNR beamforming to these sets separately to obtain the two 
beamformer outputs. 



3 Numerical Results 

The following results pertain to the APERTIF prototype system which includes a phased array feed 
consisting of 144 tapered slot antenna elements, 144 low noise amplifiers and an active (digital) beamforming 
network. The entire receiving system has been modeled by an enhanced version of the method of moments 
used as the electromagnetic simulation software, combined with a microwave circuit simulator (CAESAR), 
and the high frequency electromagnetic solver GRASP for modeling the refiector antenna O [6] . 

Figure [2l[a) illustrates the magnitude of the Jones matrix element | Jii| for the optimal method. Element 
|Jii| indicates the beam intensity that is referenced to output 1 of the Jones polarimeter. The results are 
shown for 37 beams and have been overlayed (not summed) to ease the comparison between the beams 
within the field of view. Figure [21(b) shows | Jii| along the = 0^ direction (for the on-axis and 3x2 off-axis 
beams) for the four beamforming scenarios. The values of |Jii| are the same for the optimal, eigenvalue 
and max-SLNR methods. The bi-scalar method measures a lower field intensity since half of the antenna 



output signals are combined into an x beam, and the other half into a y beam, respectively. Figure [2jc) plots 
the Jones matrix element |Ji2| for 37 beams. The element |Ji2| measures the degree of cross-polarization 
leakage for the optimal beamformer. Clearly, the polarization purity is very high in the beam centers, but 
deteriorates towards the periphery of the beam and can become as high as -15dB for the beams at the edge 
of the FoV. Figure [2]^d) shows the cut = 0^ for | J12I for the four beamforming scenarios. It is remarkable 
to see that the eigenvector beamformer leads to a cross-polarization leakage which can be as high as -5dB, 
while the cross-polarization leakage for the beams obtained with the other beamforming schemes is as low 
as -20dB for the same beam directions. Figure [2](e) indicates that for the bi-scalar method, the polarization 
homogeneity across the beam can be significantly poorer than that of the full-polarization beamformers, and 
the cross-polarization leakage can rapidly degrade to -12dB for out-of-center beams. Figure [Sff) illustrates 
the beam sensitivities at the two beamformer outputs for the optimal and bi-scalar method. In addition, the 
sensistivity bounds are shown for the optimal, bi-scalar, eigenvalue, and max-SLNR beamforming methods. 
These upper and lower bounds are the same for all methods, except for the bi-scalar method, for which 
these bounds are lower because this scheme does not beamform all antenna output signals and therefore 
compromises the sensitivity and beam polarization purity, in particular for out-of-center beams. 

Clearly, |Ji,2| ^ at the beam centers for both the optimal and max-SLNR beamformers. Hence, the 
IEEE definitions for cross-polarization, i.e. the cross-polarization discrimination (XPD) and isolation (XPI) 
values, tend to infinity at the beam centers. This surely meets the SKA specification, which states that 
the raw cross-polarization level (before calibration) must be lower than -20dB at zenith [7J. However, the 
polarimetric homogeneity of the PAF beams is not constant over the FoV and rapidly degrades as a function 
of off-axis angle of observation, especially for the bi-scalar beamformer. This implies that further calibration 
must be done to correct for the direction dependent artefacts in inter ferometric images, in particularly when 
the dominant source moves away from the centre to the periphery of the beam [8j. 
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Figure 2: Simulated Jones matrix elements for several PAF beams and polarimetric beamforming methods. 
Since the half-power-beamwidth is equal approximately to 0.5°, the beam separation distance is about one 
half- power beamwidth. 



